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We study strong solutions of the simplified Ericksen-Leslie system modeling compressible 

nematic liquid crystal flows in a domain f2 C M'^. We first prove the local existence of unique 

strong solutions provided that the initial data pQ,UQ,dQ are sufficiently regular and satisfy a 

natural compatibility condition. The initial densitiy function po niay vanish on an open subset 

(i.e., an initial vacuum may exist). We then prove a criterion for possible breakdown of such a 

Qh ' local strong solution at finite time in terms of blow up of the quantities and llVdlLs^oo. 

< .... 

'5 ■ 1 Introduction 

■ 

Nematic liquid crystals are aggregates of molecules which possess same orientational order and are 
made of elongated, rod-like molecules. The continuum theory of liquid crystals was developed by 
Ericksen [S] and Leslie [30j during the period of 1958 through 1968, see also the book by de Gennes 
^ I [II]- Since then there have been remarkable research developments in liquid crystals from both 

QQ ' theoretical and applied aspects. When the fluid containing nematic liquid crystal materials is at 

, rest, we have the well-known Ossen-Frank theory for static nematic liquid crystals, see Hardt-Lin- 

Kinderlehrer |13] on the analysis of energy minimal configurations of namatic liquid crystals. In 
■ general, the motion of fluid always takes place. The so-called Ericksen-Leslie system is a macroscopic 

continuum description of the time evolution of the materials under the influence of both the flow 
velocity field u and the macroscopic description of the microscopic orientation configurations d of 
rod-like liquid crystals. 

When the fluid is an incompressible, viscous fluid, Lin [19] first derived a simplified Ericksen- 
Leslie equation modeling liquid crystal flows in 1989. Subsequently, Lin and Liu |20l [21] made 
some important analytic studies, such as the existence of weak and strong solutions and the partial 
regularity of suitable solutions, of the simplifled Ericksen-Leslie system, under the assumption that 
the liquid crystal director field is of varying length by Leslie's terminology or variable degree of 
orientation by Ericksen's terminology. 

When the fiuid is allowed to be compressible, the Ericksen-Leslie system becomes more com- 
plicate and there seems very few analytic works available yet. We would like to mention that very 
recently, there have been both modeling study, see Morro |31) . and numerical study, see Zakharov- 
Vakulenko |39j . on the hydrodynamics of compressible nematic liquid crystals under the influence 
of temperature gradient or electromagnetic forces. 

This paper, and the companion paper [18j, aims to study the strong solutions of the flow of 
compressible nematic liquid crystals and the blow up criterions. 
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Let C M"' be a domain. We will consider the simplified version of Ericksen-Leslie system 
modeling the flow of compressible nematic liquid crystals in fi: 



pt + V ■ {pu) = 0, 
put + pu-Vu + V{P{p)) = Cu-Vd- Ad, 
dt + u-Vd = Ad+ \Vd\'^d, 



(1.1) 
(1.2) 
(1.3) 



where p : Q x [0, +oo) — t- M"*^ is the density function of the fluid, u : $7 x [0, +oo) — )• M'^ represents 
velocity field of the fluid, P = P{p) represents the pressure function, d : ilx [0, +oo) — t- S"^ represents 
the macroscopic average of the nematic liquid crystal orientation field, V- is the divergence operator 
in R^, and C denotes the Lame operator: 



where p, and A are shear viscosity and the bulk viscosity coefficients of the fluid respectively that 
satisfy the physical condition: 



We refer to the readers to consult the recent preprint [6j by Ding-Huang- Wen-Zi for the derivation 
for the system ()l.ip - ()1.3p based on energetic-variational approaches. Throughout this paper, we 
assume that 



Notice that (jl.ip is the equation of conservation of mass, (II. 2p is the equation of linear momentum, 
and (jl.Sp is the equation of angular momentum. We would like to point out that the system 
(jl.ip - (|1.3p includes several important equations as special cases: 

(i) When p is constant, the equation (II. ip reduces to the incompressibility condition of the fluid 
(V • n = 0), and the system ()l.ip - ()1.3p becomes the equation of incompressible flow of namatic 
liquid crystals provided that P is a unknown pressure function. This was previously proposed by 
Lin [19j as a simplified Ericksen-Leslie equation modeling incompressible liquid crystal fiows. 

(ii) When d is a constant vector field, the system (ll.ip -( fL2]l becomes a compressible Navier- 
Stokes equation, which is an extremely important equation to describe compressible fiuids (e.g., 
gas dynamics). It has attracted great interests among many analysts and there have been many 
important developments (see, for example. Lions [27J, Feireisl [10] and references therein). 

(iii) When both p and d are constants, the system (|l.ip - (ll.2p becomes the incompressible Naiver- 
Stokes equation provided that P is a unknown pressure function, the fundamental equation to 
describe Newtonian fiuids (see, for example. Lions |26) and Temam [M] for survey of important 
developments) . 

(iv) When p is constant and u = 0, the system (jl.ip - (|1.3p reduces to the equation for heat fiow 
of harmonic maps into S"^ . There have been extensive studies on the heat flow of harmonic maps in 
the past few decades (see, for example, the monograph by Lin- Wang [24J and references therein). 

From the viewpoint of partial differential equations, the system (ll.ip - (jl.3p is a highly nonlinear 
system coupling between hyperbolic equations and parabolic equations. It is very challenging to 
understand and analyze such a system, especially when the density function p may vanish or the 
fluid takes vacuum states. 

In this paper, we will consider the following initial condition: 



Cu = pAu + (p + A)Vdiv u, 



p>0, 2/i + 3A>0. 



(1.4) 



P : [0, +oo) — )■ R is a locally Lipschitz continuous function. 



(1.5) 



{p,u,d) 



t=o 



{po,uo,do) 



(1.6) 



and one of the three types of boundary conditions: 



2 



(1) Cauchy problem: 

= M^, and p, u vanish at infinity and d is constant at infinity (in some weak sense). (1-7) 

(2) Dirichlet and Neumann boundary condition for (u, d): 17 C is a bounded smooth domain, 



dd, 
ov 



and 

rid , ^ 

0, 1.8 

y an 

where v is the unit outer normal vector of dVl. 

(3) Navier-slip and Neumann boundary condition for (u,d): C is bounded, simply con- 
nected, smooth domain, and 

dd 

(u • I/, (V X n) X ly, — ) = 0, (1.9) 

Oh' dU 

where V x u denotes the vorticity field of the fluid. 

To state the definition of strong solutions to the initial and boundary value problem (ll.l[ )-( fL3]l . 
(|1.6p together with (|1.7p or (jl.Sp or (jl.9p . we introduce some notations. 

We denote 

f dx = / / dx. 



For 1 < r < oo, denote the spaces and the standard Sobolev spaces as follows: 

U' = D^'"- = {-u G L\^^{^) : WV^uWir < oo} , 

L)i = jii G : ||V'u||i2 < oo, and satisfies ([TTl) or ([THl) or ([LS]) for the part of n|, 

Hl = L^r\Dl, ||n||^fe,. = IIV^'uIIl'-. 



Denote 
and let 



Qt = VLx [0,r] (T>0), 



denote the deformation tensor, which is the symmetric part of the velocity gradient. 

Definition 1.1 For T > 0, {p,u,d) is called a strong solution to the compressible nematic liquid 
crystal flow M.l\j - [L3^) infix {0,T], if for some q G (3,6], 

< p G C([0,r];VFi'^n^'), Pt G C([0,T];L2nL''); 

Vde C{[0,T];H^)f]L\0,T;H^), dt e C{[0,T]; H^)f] L^{0,T; H^), \d\ = 1 in Q^; 
and {p,u,d) satisfies il.l\ }- [T73\} a.e. inVtx (0,T]. 

The first main result is concerned with local existence of strong solutions. 
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Theorem 1.2 Assume that P satisfies iL^, po > 0, po £ W^'" C]H'^ C]L'^ for some q E (3,6], 
uq G D^f]DQ, Vdo € and \do\ = 1 inQ. If, in additions, the following compatibility condition 

Cuo - V{P{po)) - Ado ■ Vdo = y/pog for some g G L^{n, M^) (1.10) 

holds, then there exist a positive time Tq > and a unique strong solution {p,u,d) of il.l\) - [T73\) . 
(EJj together with {T^j or or [L^ in Q x {0,To]. 

We would like to point out that an analogous existence theorem of local strong solutions to the 
isentropic0 compressible Naiver-Stokes equation, under the first two boundary conditions (jl.7p and 
(II. Sp . has been previously established by Choe-Kim [4j and Cho-Choe-Kim [3]. A byproduct of 
our theorem 11.21 also yields the existence of local strong solutions to a larger class of compressible 
Navier-Stokes equations under the Navier-slip boundary condition (jl.9p . which seems not available 
in the literature. 

In dimension one, Ding-Lin- Wang- Wen |7j have proven that the local strong solution to (II. ip - 
(jl.Sp under (|1.6p and (jl.Sp is global. For dimensions at least two, it is reasonable to believe that 
the local strong solution to (|l.ip - (|1.3p may cease to exist globally. In fact, there exist finite time 
singularities of the (transported) heat flow of harmonic maps ()1.3p in dimensions two or higher (we 
refer the interested readers to [23| for the exact references). An important question to ask would 
be what is the main mechanism of possible break down of local strong (or smooth) solutions. 

Such a question has been studied for the incompressible Euler equation or the Navier-Stokes 
equation by Beale-Kato-Majda in their poineering work pLj, which showed that the L°°-bound of 
vorticity V x u must blow up. Later, Ponce [29] rephrased the BKM-criterion in terms of the 
deformation tensor T>{u). 

When dealing with the isentropic compressible Navier-Stokes equation, there have recently 
been several very interesting works on the blow up criterion. For example, if < T^, < +oo is the 
maximum time for strong solution, then (i) Huang-Li-Xin [15] established a Serrin type criterion: 
limrtr. {\\di^u\\L^o,T;L°-) + \\Vpu\\l''{o,T;L-)) = oo for | + ^ < 1, 3 < r < oo; (ii) Sun- Wang-Zhang 
[35], and independently [T5], showed that if 7p > A, then limT^T^, IIpIIl°°(o,T;L°°) = oo; and (iii) 
Huang-Li-Xin [16j showed that liniT^T* II^(^)IIli(o,T;L°°) = co- 

When dealing the heat flow of harmonic maps (|1.3p (with n = 0), Wang [5B] obtained a Serrin 
type regularity theorem, which implies that if < < +oo is the first singular time for local 
smooth solutions, then lim^-^j'^^ ll^^llL2(o,r;L°°) = c». 

When dealing with the incompressible nematic liquid crystal flow, Lin-Lin- Wang [25j and Lin- 
Wang [23] have established the global existence of a unique " almost strong" solutioro for the initial- 
boundary value problem in bounded domains in dimension two, see also Hong [14] and Xu-Zhang 
[38] for some related works. In dimension three, for the incompressible nematic liquid crystal flow 
Huang- Wang [17] have obtained a BKM type blow-up criterion very recently, while the existence 
of global weak solutions still remains to be a largely open question. 

Motivated by these works on the blow up criterion of local strong solutions to the Navier-Stokes 
equation and the incompressible nematic liquid crystal flow, we will establish in this paper the 
following blow-up criterion of breakdown of local strong solutions under the boundary condition 
(ini) or (1121). 

Theorem 1.3 Let {p,u,d) be a strong solution of the initial boundary problem ^.1\ )- [T7^) . ^.6\) 
together with (fi. 7\ ) or il.8\) . Assume that P satisfies il.5]) . and the initial data {po,uo,do) satisfies 



i.e. P = ap'' for some o > and 7 > 1. 

that has at most finitely many possible singular time. 



4 



U.10\) . If < < +00 is the maximum time of existence and 7/i > 9A, then 



lim 



( 



P\\l°°{q, 



>,T;L°°) + l|V(i||L3(o,T;L°°) 



= 00. 



(1.11) 



We would like to make a few comments of Theorem 11.31 

Remark 1.4 (a) Since we can't yet prove Lemma 4.2 for the Navier-slip and Neumann boundary 
condition (|1.9p . it is unclear whether Theorem 11.31 remains to be true under the boundary condition 



(b) In [18j . we obtained a blow-up criterion of (jl.ip - (jl.3p under the initial condition (jl.6p and the 
boundary condition (II. 7p or (II. Sp or (II. 9p in terms of u and Vd: if < T^, < +00 is the maximum 
time of existence of strong solutions, then 



(b) For compressible liquid crystal flows without the nematicity constraint {\d\ = ljf|, Liu-Liu |22] 
have recently obtained a Serrin type criterion on the blow-up of strong solutions. 

(c) It is a very interesting question to ask whether there exists a global weak solution to the 
initial-boundary value problem of (ll.ip - (ll.3p in dimensions at least two. In dimension one, such an 
existence has been obtained by Ding- Wang- Wen 

Now we briefly outline the main ideas of the proof, some of which are inspired by earlier works 
on the isentropic compressible Navier-Stokes equations by [3], [35], and [T6]. To obtain the existence 
of a unique local strong solution to (11.11 - (|1.3p . under (|1.6p and (]1.7p or (jl.Sp or (]1.9p . we employ 
the Galerkin's method that requires us to establish a priori estimate of the quantity 



for strong solutions (p, n, d) in the form of a Gronwall type inequality. See Theorem 2.1. It may be 
of independent interest that we establish IF^''^-estimate for the Lame equation under the Navier-slip 
boundary condition, see Lemma 3.1. 

To prove the blow-up criterion (jl.lip of Theorem 11.31 in terms of p and Vd, a critical step is 
to establish the L^L^-estimate of Vp. From the continuity equation (II. ip . this requires that the 
Lipschitz norm of velocity field u, or ||V^ii(t)||i9 is bounded in Lj. This is done in several steps. 

(1) We show that under the condition Jfi > 9A, the bound of (||/o||L^Lg° + II ™ equations 
(II. 2p and (II. 3p can yield both a high integrability and a high order estimate of u and Vd, i.e. both 
(II/'^'"IIlj°°l5 + ||V(i||j;^oo^5 ) and (||Vu||2,^oo^2 + || V^(i||j;^oo^2 ) are bounded. See Lemma 4.2. 

(2) Based on these estimates from (1), we establish that V^d is bounded in Lf°L'^ and Vn is 
bounded in L^W^''^ + Lf^(BMOi.). To achieve it, we adapt the approach, due to Sun- Wang-Zhang 
[35] . by decomposing u = w + v, where v G Hq{^) solves the Lame equation Cv = V{P{p)). 
One can prove that Vv £ L^(BMOa;) by the elliptic regularity theory. The difficult part is to 
show that V'^w G L^L% for 3 < q < 6. In order to obtain this estimate, we first establish that 
(||->/p'u||j;^ooj;^2 + 11 Vdt ||iooj;^2 ) and (||Vtt||j;^2j^2 + ||(i4t||2,22,2 ) are bounded by viewing (|1.2p as an evolution 
equation of the material derivative ut + u ■ Vu and performing second order energy estimates 
of both equations (II. 2p and (jl.3p . Then we employ VF^''^-estimate of the Lame equation to control 
IIV^wIIl?. The details are illustrated by Lemma 4.4 and Corollary 4.5. 

^the right hand side of equation (|1.3|l is replaced by Ad + f{d) for some smooth function / : R'^ — >■ R'^, e.g. 



dLSD- 




Hin M|P('u)||ii(o,r;L°°) + l|V(i||i2(o,T;L-') ) =+00. 



Mt)\\mnw^-^ + l|V^x(t)||i,2 + \\^ut{t)\\L2 + \\V^d{t)\\L2, 3 < g < 6 



/(d) = (|d|2 - l)d. 
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(3) We show that ||V/?||j;^2pi2,'j is bounded by an argument similar to [35] §5. Then we apply 
VF^'''-estimate of the Lame equation again to control || V^u||2,^2,2 and ||n||^oo£,2,<7. See Lemma 4.6, 
Corollary 4.7, and Corollary 4.8. 

It is interesting to notice that during the proof of both the existence of a unique local strong 
solutions and the blow-up criterion for strong solutions, specific forms of the pressure function P{p) 
play no roles and it is the local Lipschitz regularity of P that matters. 

The paper is written as follows. In §2, we derive some a priori estimates for strong solutions or 
approximate solutions via the Galerkin's method. In §3, we prove both the local existence by the 
Gakerlin's method and uniqueness of strong solutions. In §4, we discuss the blow up criterion of 
strong solutions and prove Theorem 11.31 

Acknowledgement. The first two authors are partially supported by NSF grant 1000115. The 
work is completed during the visit of third author to University of Kentucky, which is partially 
supported by the second author's NSF grant 0600162. The third author wishes to thank the 
department of Mathematics for its hospitality. 

2 A priori estimates 

In the section, we will derive some a priori estimates for strong or smooth solutions {p, u, d) to (jl.ip - 
(jl.3p on a bounded domain, associated with the initial condition (|1.6p and the boundary condition 
(jl.Sp or (jl.9p . provided that the initial density function has a positive lower bound, po ^ > 0- 
All these a priori estimates we will obtain are independent of 6 > and the size of the domain 
when Q = Br (R > 1) is a ball in M^, which are the crucial ingredients to prove the local existence 
of strong solutions to (ll.ip - (11.3p when we allow the initial data po > and unbounded domain 
$7 = M^. Although these estimates may have their own interests, we mainly apply them to the 
approximate solutions to (jl.ip - (jl.3p that are constructed by the Galerkin's method. 

Throughout the paper, we denote by C generic constants that depend on UpollvFi-'Jnii'inLi > 
ll^ollD^nDi' ll^^^olliTa, and P, but are independent of 5 > 0, the solutions {p,u,d) and the size of 
domain when Q = Bji {R> 1) is a ball in M^. We will also use the obvious notation 

k 

II • lUin-nXfe = X] II ■ 11^' 
1=1 

for Banach spaces Xi, I < i < k and A; = 2,3. We will use A B to denote A < CB for some 
constant generic C > 0. 

Let u, d) be a strong solution of (|l.ip - (|1.3p in 17 x (0, T] (or the approximate solutions 
{p'^,u"^,d^) of ()l.ip - (jl.3p constructed by the Galerkin's method in §3.2 below). For simplicity, we 
assume < T < 1. For < t < T, set 

:= sup (\\p{s)\\Hinw^,. + l|Vn(s)||i2 + Mut{s)\\L2 + \\V^d{s)\\Hi + l) • (2.1) 

0<s<t ^ ' 

The main aim of this section is to estimate each term of $ in terms of some integrals of In §3 
below, we will apply arguments of Gronwall's type to prove that $ is locally bounded. 

Throughout this section and §3, we will let T to denote the set that consists of monotonic 
increasing, locally bounded functions M from [0, +oo) to [0, +oo) with M(0) = 0, which are inde- 
pendent of b and the size of 0. The reader will see that the exact form of M € -F is not important 
and may vary from lines to lines during the proof of the Lemmas. 

Now we state the main theorem of this section. 
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Theorem 2.1 There exists M £ T such that for any < t < T , it holds 



where 



$(t)<exp CM{po,uo,do) + C / A/(^>(s))ds 
L Jo 

£uo-V(P(po))-Ado-Vdo 



M{po,uo,do) = 1 + 



(2.2) 



(2.3) 



L2 



The proof of Theorem 12.11 is based on several Lemmas. We may assume P(0) = 0. Observe 
that (II. 5p imphes that the Lipschitz norm 



Bp{R) := ||-P'||l=o([o,_r]) : [0, +cxd) — )■ [0, +oo) is montonic increasing and locally bounded. (2.4) 
Lemma 2.2 (energy inequality) There exists M G such that for any < t < T , it holds 



Jn 



+ |V(i|2) dx+ / / iVul"^ + \Ad+\Vd\'^d\^ dx<C+ I M($(s))ds. (2.5) 



Proof. Here we only sketch the proof for the boundary condition (11. 9p . Multiplying (11. 2p by u 
and integrating over Q, using An = Vdivn — V x (V x ti) and (jl.ip . and applying integration by 
parts several times, we obtain 

^^yp|'up(ix + y"(^|V X lip + (2^ + A)|div'up)da; = J P{p)dwudx - J u -Vd ■ Addx. (2.6) 

Since Q is assumed to be simply connected for the boundary condition (|1.9p . we have (see j37]): 

I|Vu||l2 ^ l|V X u\\l2 + ||divti||i2, Vn G i7^(Q) with n • = on Sfi. (2.7) 
This and (jl.ip imply 

y (^|V X -up + (2^ + A)|divn|2) c^x > | y(|V x up + jdivup) dx > ^ J jVup dx. (2.8) 



By Cauchy inequality, we have 



P(p)divn dx 



-2^/ + ^ / \Pip)\'dx. 



(2.9) 



Multiplying (|1.3p by Ad + |Vdpd and integrating over Q, using integration by parts and the fact 
that \d\ = 1 we obtain 



ld_ 
2di 



j iVdp dx + j \Ad+ iVdpdp dx = j u-Vd- Addx. 
Combining i^M), <^M, '^M, and (IXTUD together, we obtain 

^ J\p\u\'^ + \Vd\'^)dx + Ji^lVul"^ + \Ad+\Vd\'^d\'^)dx <C J \P{p)\'^ dx. 

To estimate the right hand side of (j2.1ip . first observe that by (|2.4p we hav^ 

IIpIIl- + \\Pip)\\L^ + \\Pip)\\mnw^'^ <C^ + CBpi\\p\\L^)<i> < M($) 

*when Q, — Bb. for i? > 1, one can the independence of C with respect to R as follows: 

\\p\\l-^{Br) < max ||p||l~(Si(x)) < C max \\p\\w^,^Bi{x)) < C\\p\\w^.i(Bii)- 



(2.10) 



(2.11) 



(2.12) 
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for some M €z J-. It follows from (jl.ip and Sobolev's inequality that 

J \P{p)\^dx = J \P{po)\^dx + 2 J P{p)P'{p){-pdwu-Vp-u)dxdt 

< C + C ri?p(||HU^)(||P(p)||i3||VH|i2 + ||P(p)||,.2||/5|U^)||Vn||i2ds 

Jo 

< C+ f M{^{s))ds <C + M{^{t)) (2.13) 

JO 

as M(<I>(s)) is increasing and t <1. Substituting (|2.13p into (j2.1ip and integrating over [0,t] yields 
(I23D . □ 

Now we want to estimate ||Vn(t)||^i in terms of 
Lemma 2.3 There exists M ^ T such that for < t < T, it holds 

||Vn(t)||^,i <Mm)). (2.14) 

Proof. By the standard i/^-estimate of the Lame equation with respect to the boundary condition 
(fLT]) or ^Bj or ([131), (1212]), and Holder's inequality, we have 



\Vu\\jJ^ <\\^U\\12 + \\V 



u\ 



L2 



<\\put\\l2 + \\pu ■ Vu\\l2 + \\y{P{pmh + 11^^^ • Vd||i2 + llV^xll^^ 

<\\p\\l^ WVpMI^ + ll/'lli- \M% llV^^llia + S|>(||/9||L-)||Vp||i2 (2.15) 

2 ||V7JI|2 , Ilv7„.ll2 



+ \\Ad\\i4Vd\\ie + ||Vn||£2 
<M($)(1 + ll^lliellVi^llia) + C\\Ad\\UVd\\l, 

for some M € J^. By the interpolation inequality, Sobolev's inequalitjil, we obt, 



am 



II«IIl6||V"IIl3 <C||Vn||'^2||Vu||^^i. (2.16) 
Similar to (j2.16p . by (|2.5p . we obtain 



< 



\\Ad\\UVd\\l,<\\AdU4M\L4^d\\l. 
\\Ad\\U\^d\\l, + \\Ad\\l4VM\l,<M{<P) 



(2.17) 



for some M £ T. Substituting (I2.16p . (|2.17p into (|2.15p . and using (j2.5p and Cauchy's inequality, 
we have 

||Vtx||^i <^||Vu||^i +M{^{t)) 
for some M £ J^. This gives (j2.14p and completes the proof. □ 
Now we want to estimate ||-y//0'Ut||L2- More precisely, we have 



^ when Q = Bjt for > 1, by simple scalings, one has 
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Lemma 2.4 There exists M G J- such tha for any < t < T , it holds 



\Vut\^dxds < CM{po,uo,do) + [ M{^{s))ds. 

Jo 

Proof. Differentiating (|1.2p with respect to t, we 



p\ut\ dx - 
n Jo Jn 



(2.18) 



putt + pu ■ Vut + ptut + ptu ■ Vu + put ■ Vn + V{P{p))t 

=(2^ + A)Vdivuj - pV X (y X ut) -V ■ (Vdt Vd + Vd(g) Vdj - Vd-Vdt h). 



(2.19) 



Multiplying (|2.19p by Ut, integrating the resulting equations over Q, and using (jl.ip and integration 
by parts, we have 



ld_ 
2di 



J p\ut\'^dx + j ((2/i + A)|divutp + /i|V X -Ufp) dx 
= — 2 j puut ■ Vut dx — j ptu ■ Vn ■ utdx — j put ■ Vu ■ utdx + J P'{p)ptdiv ut dx 

5 

+ / {Vdt Vd + Vd Vdt - Vd ■ Vdt I3) : Vut dx = Ih. 

i=l 

By Holder's inequality, Sobolev's inequality, (I2.12p . and (I2.14p . we have 

l-^^ll S\^Ut\\LA\\fpUt\\LA\\fpAL^ 

<\Wut\\L2\\^Ut\\LA\Vp\\L^\\^'^\\m < M{^)\\Vut\\L2 

for some M £ J^. 

By (jl.ip . Holder's inequality, Sobolev's inequality, (I2.12I) . and (|2.14p . we have 



\Ih 



J pu ■ V(u • Vu ■ Ut) dx 



= J P'^ ' ■ • ut + u ■ VVu ■ ut + u ■ Vu ■ Vut) dx 

+ ll\/^l|L°°lkllL6||V'u||^6||VP^illL2 

<M{^){l + \\Vut\\L2) 



for some M £ J^. For I/3, by (j2.14p we have 

l^-^sl <ll\/p||L°°|IVpMt||i2||Vu||i3||nt||i6 

<ll^/p||L°o||Vpni||i2||V^z|||,||Vn|||,||Vn^||i2 <M($)||Vni||i2 
for some M G J". For Ih, by ([Ll]), (|242]) . and (f23ip we have 

1^-^41 <Bp{\\p\\L^)\\pt\\L2\\diY UtWi^ 

<Bp{\\p\\loo){\\V p\\l2\\u\\l°° + ||p||Lo°||divii||^2)||divMt||j^2 
<Mm\Vuth2 



%eie we have used the fact that Ad- Vd = V ■ (Vd® Vd- ilVdj^Ia), where Vd® Vd = {d^^ ■ dx^)_^ 
is the identity matrix of order 3. 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



<i,j<3 



and I3 
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for some M E 7^. For I/5, by (j2.5p we have 



l^^5|< / \Vd\\Vdt\\Vut\dx<\\Vut\\L4'^d\\L^\\Vdt\\L2 

<\\VUt\\L2\\Vd\\H2\\Vdt\\L2 ^^-^^^ 
<\\'^Ut\\L2i\\Vd\\L2 + \\VmHi)\\Vdt\\L2 < {C + Mm\Nut\\L4^dt\\L2 

for some M £ T . Substituting (j2.2ip - (j2.25p into (j2.20p . and using Cauchy's inequality, we have 

<^ / |V«,l^<ii + Af(*) + (C + Af(*))||Vd,|||, 

for some M G J^, where we have used the following inequality due to ^37j: if (i) either is simply 
connected and w • = on 90 or (ii) n = on dVl |j, then 

||Vut||i2 < ||divnt||i2 + ||V X utWiji. {2.21) 

By (|2:26]) . we have 

^ y" plntl^dx + i y" iVu^pdx < M($) + (C + M{^))\\Vdt\\l2. (2.28) 

Differentiating (jl.3p with respect to x, we have 

Vdt - VAd = V(|Vdpd) - V(u • Vd). (2.29) 

From ([2:29]) . we have| 

l|Vdt|U2 <||Vn • Vd|U2 + lit. • V^d\\L2 + ||VAd|U2 + llVdIlie + ||Vd • V^dl^a 
<||Vd||i^ ||V^/||i2 + 11^x11^6 llV^dll^a + + (1 + ||V2d||^2)3 

+ ||V(i||L-||V2d||i2 

<||Vd||H2||Vn||i2 + \\Vu\\L2\\V^d\\Hi + ||VAd||i2 + (1 + \\V^d\\L2f 

+ \\Vd\\H2\\V^d\\L2 

<M($) + 1 

for some M £ J^. 

Substituting (I2.30p into (I2.28j) . and using Cauchy's inequality, we have 



d_ 
lit 



J p\ut\'^dx + J J \Vut\^dxds < C J p\ut\^ dx ^ + j M{^{s))ds + C 



rt 



< CM{po,uo,do)+ / Mms))ds 
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(2.30) 



■ j p\ut\'^dx + ^ j \Vut\'^dx < M($) + C (2.31) 
for some M G J^. Integrating (j2.3ip over (0,t), and using (|1.2p . and (jl.lOp . we have 



^in fact, in this case, the inequahty (|2.27p is an equahty. 

*here we also use the Sobolev's inequality: ||Vd||£oo(Q) < C[| Vci||£f2(f2) and the fact that C can be chosen inde- 
pendent of 7? when Q = Br for R>1. 
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for some M €z J-. This completes the proof. □ 

As an immediate consequence of Lemma 12.41 we obtain an estimate of ||Vn||j;^2. 
Lemma 2.5 There exists M ^ T such that for < t < T , it holds 

J \Vu{t)f dx < CM{po,uo,do) + M{^{s))ds. (2.32) 

Proof. By Cauchy's inequahty, Lemma l2.2p . Lemma 12.31 and Lemma 12.41 we have 

f \Vu\'^{t)dx= I \Vuo\'^dx + 2 [ [ Vu-Vutdxds 
J J Jo Jn 



<C + I I \Vu\^dxds+ I I \Vut\^dxds 
Jo Jn Jo Jn 

<CM{po,uo,do)+ I M{^is))ds 
Jo 



for some M £ This completes the proof. □ 
Lemma 2.6 There exists M £ such that for < t < T, it holds 

do) + C Mms))dsY (2.33) 

Proof. It follows from [3] (page 249, (2.11)) that 

\\p{t)\\H^nw^-i <\\po\\mnw^-ieW^C WVuWHinD^.ids'^ ■ (2-34) 

By W^^'^-estimate of the Lame equation under either Dirichlet boundary condition (jl.Sp or the 
Navier-slip boundary condition p.9|) (see Lemma l3.ll below). (|1.2|) . and Sobolev's inequality, we 
have 

4 

IIV^uIIl, <\\pMli + \\pu ■ Vu\\l, + ||V(P(p))||l, + \\Vd ■ Ad^, = ^^^i- (2-35) 

i=l 

If g = 6, then by Sobolev's inequality we have 

Ilh < ||p||L-||nt||i6 <$||Vnt||i2. (2.36) 
If g G (3, 6), then by Holder's inequality and Sobolev's inequality, we have 

nil < M^^Wnth^ < IIHiyi|p|lL^^I|Vnt||i2 <$||Vnt||^2, (2.37) 

where we have used the fact that J pdx = J podx. From ()2.36p and ()2.37p . we have that for 
(3,6], 

IIh<<f\\Vuth2. (2.38) 



11 



For III2, if ? G (3, 6], then by similar arguments, Lemma 2.2, and Lemma 2.3, we have 

III2 < ^||Vn||^i < M($) (2.39) 
for some M € J^. For III3 and III4,, if g E (3, 6], then we have 

Ilh + Ilh < CBpi\\p\\L^)\\Vp\\L. + ||Vd||^2 < M($) (2.40) 
for some M £ Substituting ([235]), and H^IMh into ([2321), we have 

llVnlli, <nVut\\L2 + M($) < llVn^lli^ + M($) (2.41) 
for some M £ F. Integrating (j2.4ip over (0,i), and using Cauchy's inequahty and (j2.18p . we have 

f ||V\||l, < CA^(po, uo, do) + r M($(s)) ds. (2.42) 
Substituting (|2TH) and (I2:i2]l into ([23i|) . we have 

\\p{t)\\mnw^^^ <expl^CM{po,uo,do) + C M($(s))ds| 

for some M €z J-. This completes the proof. □ 
Lemma 2.7 There exists M £ F such that for any < t < T , it holds 

\N^d\\h+ [ \Ndt\\l2ds<C+ [ M{^{s))ds. (2.43) 
Jo Jo 

Proof. Multiplying (j2.29p by Vdf and integrating over using integration by parts and ^ = 
on d^l, we obtain 

\\Vdt\\l2 + ~\\Ad\\l2= J [V{\Vd\^d)-V{u-Vd)] Vdtdx 

<\\\'^dt\\l2+C j |V(|Vdpd)pdx + C7 j \V{u-Vd)\^ dx. 

Thus we have 

||Vdt||i. + ^||Ad||i, < j |V(|Vdpd)|2dx + j \V{u-Vd)\''dx. (2.44) 
Similar to the proof of (|2.30p , we obtain 

||Vdt||i, + |||Ad||i, <M($) (2.45) 



for some M £ T. Integrating (j2.45p over (0,t) and applying W^^'^-estimate of the equation (jl.3p . 
we have 

l|V^ti|li2+ / llVdtlliads < ||V2(io||i2 + / M($(s))ds <C7+ I M{<^{s))ds. 
Jo Jo Jo 

This completes the proof. □ 
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Lemma 2.8 There exists M G T such that for < t < T, it holds 

||v3d||2, +^*||V2dt||2,ds < (^(:7A^(po,uo,4) + ^*M($(s))ds^ . (2.46) 

Proof. Multiplying (j2.29p by VAd^, integrating over 0,, using ^ = on dO, and integration by 
parts, we obtain 

-IT p 

\\^dt\\l2 + --f\\VAd\\l2 = / [V{u-Vd) - V{\Vd\^d)] - VAdtdx 

=^ / [V(u- Vd) - VdVdpd)] • VAddx (2.47) 

~ \ it ^^"^ ' ^'^^ ~ ^(l^^l^^)] ■ VAddx. 
Now we need to estimate the second term of right side as follows. 

- J -^[yiu-Vd)] ■VAddx = - J [Vut ■ Vd + Vu ■ Vdt + ut ■ V'^d + u ■ V'^dt] -VAd 

1=1 

By Holder's inequality and Sobolev's inequality, we have 

\IVi\ < \\Vut\\L4'^d\\L^\\VAd\\L-2 < \\Vut\\L2\\Vd\\j,2<M{^) + \\Vut\\h (2.49) 
for some M €z J-. 

By Holder's inequality, Sobolev's inequality, (|2.14|) . (|2.3U|) and Young's inequality, we obtain 



dx 

(2.48) 



(2.50) 



(2.51) 



im <||Vu||i6||Vdt|U3||VAd|U2 

<||Vu||j^i||Vdt|ki||VAd|U2 

<||Vn||j^i||v2dt||i2||VA(i||i2 + || Vn||j^i ||V(it [1^2 ||VA(i||i2 
<e\\V^dt\\l2 + M{^) 

for some M £ T . 

By Holder's inequality, Sobolev's inequality and Cauchy's inequality, we obtain 

l^^sl <||t^t||L«l|V2d||i3||VAd||i2 

<||V^Xi|U2||V2d||^i ||VAd||i2 < M{^) + ||Vnt||i2 

for some M £ T . 

By Holder's inequality, Sobolev's inequality, (|2.14p and Cauchy's inequality, we obtain 

iml <||n||i^||v2dt||i2||VAd||i2 < ||V7x||^i||v2dt||^2||VAd||i2 
<e||V2(it||22 +M($) 

for some M £ J-. 

Combining (H^H]), (E^H), ([230]) . ([231]) and (l232]l . we obtain 

— [V(u • Vd)] ■ VAddx <2e\\V^dt\\l2 + C\\Vut\\l2 + M($) (2.53) 



(2.52) 
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for some M E T. 

By Leibniz's rule and the fact \d\ = 1, we have 

^[VdVdl^d)] . VAddx 



< 



[\Vd\^\Vdt\ + \Vdt\\V'^d\ + \Vd\\V^dt\ + \V d\\V'^ d\\dt\]\V Ad\ dx = ^Vi. 



(2.54) 



i=l 



By Holder's inequality, Sobolev's inequality and (I2.30p . Cauchy inequality, and Young inequality, 
we obtain 

l^il <||Vd||ioo||Vdt||i2||VA(i|U2 < ||Vd||^2||Vdt||i2||VAd||i2 < M{^), (2.55) 

72^11 IIV7AJII <^ IIV7J II IIV72j||2 



\V2\ <\\^dthe\\V'dhs\\VAdh2 < \\Vdt\\Hi\\V'drm 

72j II . I IIV7J II ^^ ^ ^IIV72j ||2 



(2.56) 
(2.57) 
(2.58) 



(2.59) 



<^{\\V'dt\\L2 + ||Vdt||i2) < e\\V'dt\\i2+M{'^), 
\V3\ <||Vd|U^||V2dt||i2||VAd||i2 

<||Vd||H2||v2di||^2||VAd||^2 < e\\V^dt\\l, + M{'^), 

\Va\ <\\dt\\L4^d\\L^\\V^dh4^Ad\\L2 

<\\dt\\HANd\\H2\\V^d\\H4VAdU2 < \\dt\\H^Mi^) 
for some M €z J-. Notice that 

\\dth2<\\Adh2 + \\Vd\\l, + \\u-Vdh2 

<l|Vd||^i + ||n||i6||Vd|U3 + 1 < llVdll^i + \\Vu\\L2\\Vd\\H. + 1 < 

Thus by ([230]) . (1238]) and ([239]) . we have 

|^4| < M($) (2.60) 
for some M £ T. Combining (l23iD . (I235D . (I236D . (12371) and (l230|) . we have 

/ ^[^^l^'^l^'^)] ■ '^^^^^ ^ 2e||V2(it||22 + M(^>) (2.61) 

for some M G J^. 

Putting ([233]) and (f23T]) into (ITiTl) . we obtain 



\Adt\\l2 + -^\\VAd\\l2 <^ [ [V(n- Vd) - VdVdpd)] • VAddx 

+4e\\V^dt\\l2 + CllVntll^a + M(^>) 



(2.62) 



for some M ^ Integrating (|2.62p over (0, t), using {k = 2, 3) estimate of the elliptic equations, 
and choosing e small enough, we have 

l|V^d|li2 + [' \\V^dt\\l2ds 
Jo 

< j |V(n • Vd) - V(|Vdpd)| |VAd| dx + j |V(no • Vdo) - V(|Vdopfio)| |VAdo| dx (2.63) 

+ l|V'do|li2+ f\\Vut\\l2ds+ f M{^s))ds. 
Jo Jo 
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For the first term of right side of (j2.63p . we have 

j |V(u • Vd) - V(|Vdpd)| \VAd\ dx 



< 



(iV-ullVdl + \u\\V^d\ + \Vd\^ + iVdllV^dl) \VAd\ dx = "^VIi. 



(2.64) 



1=1 



By Holder's inequality, Nirenberg's interpolation inequality, (j2.5p . and Young's inequality, we obtain 

\Vh\ <\\^d\\L^\\Vu\\L2\\VAd\\L2 < \\Vd\\l\\Vd\\l,\\Vuh2\\VAd\\L2 

<ll Vd|||, \\Vuh. WV^Il. + II V^dlll, II Vnlli. (2.65) 
<e||V3d||i.+C(||Vd|||,||V.x||i. + ||Vn||i.), 

|V^^2| <ML4^'dh4VAdh2 < ||V^x||^2||V2d|||,||V2(i|||,||V3d||^. 

<||Vnb2||v2d||^.||V3d||,^2 + ||Vn||i.||V^d||il|V3d|||, (2-66) 
<e||V3d||i. + C||V2d||i.(||Vn||i. + ||Vn||i.), 

I^^sl <||Vd||i6||VAd||i. < ||Vd||f,i||v3d|U2 <e||v3d||i.+C7||Vd||f,i, (2.67) 



and 



l^/4| <||Vd||L^||V^(i||i2||VAd||^2 < ||Vd||2,||Vd||^,||V^d|U2||V^d|U2 

< II Vd|||. II V^dll^. II V^dll^. + II V^dlll, II V^dll^. (2.68) 



<e||V3d||i. + C(|| Vd|||, II v^diii. + iiv^diii.) 

Combining (ITeiD . (fTeSD . I^Mh . (|2:67D and I^Mh . we obtain 



j |V(n . Vd) - V(|Vd|2(i)| \VAd\ dx 



<4£||V^d||i. +C||Vd|||,,(||Vn||i. + llV^dlli.) + t^l|V^d||i.(||Vu||i. + \\Vu\\h) 
+C{\\Vd\\%, + \\vm% + \\Vu\\l,) 



(2.69) 



<Ae\\V^d\\l2 + (^CM{po,uo,do) + M{^{s))ds^ 

for some M €z J-, where we have used Lemma 2.2, Lemma 2.5, and Lemma 2.7 in the last step. 

Substituting (j2.69p into (|2.63p . choosing e small enough, and using (j2.18p . Cauchy's inequality, 
Lemma 12.51 and (j2.43p . we have 

l|V='d||i2+ f \\V^dt\\l2ds < (cM{po,uo,do)+ f Mms))ds 
Jo \ Jo / 

for some M £ J^. This completes the proof. □ 

Proof of Theorem 12.11 It is readily seen that the conclusion follows from (j2.18p . (j2.32p . (j2.33p . 
([2:43]) and (lOHD . □ 
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3 Proof of Theorem [13] 



3.1 W^'P-estimate 



In this subsection, we give a proof of VK^'^-estimate of the Lame equation on a simply connected, 
bounded, smooth domain with the Navier-shp boundary condition, which is needed in our proof of 
Theorem 11.21 We beheve that such an estimate may have its own interest. 

Lemma 3.1 For any simply connected, smooth bounded domain Q C M^, 1 < p < +00, and 
f G LP{n,R^), IfueH'^n H'^{n,R^) is a weak solution of 



Cu = / in ri, 
u-v = {yxu)^i> = on dVt. 

Then u € VF^'^(O), and there exists C > depending on p, il, and C such that 



< C 



LP 



LP + ||Vn||x,2 



(3.1) 



(3.2) 



Proof. By the duahty argument, we may assume 1 < p < 2. Since u • z/ = on di}, it follows 
from Bourguignon-Brezis that 

llV^nlliP < ||V(div u)\\lp + ||V(curl u)\\lp + \\Vu\\lp. (3.3) 

Also, since il. is simply connected and (V x n) x 1/ = on dQ, it follows from Wahl |37] that 

||V(curl u)\\lp < C\\V X curl u\\lp + ||V • (curl u)\\lp = C\\V x (curl u)\\lp 

< -\\\Cu\\LP + {2fX + X)\\V{divu)\\LP 

fi L J 

< ||V(divu)||LP + ||/||LP. 

Now we estimate ||V(div u)\\lp by the duality argument: for p' = 

||V(divn)||LP < Csupj / V{divu)-gdx : g £ C'^{n,M.^), 



(3.4) 



I^IIlp' = i}- 



For any g £ C"^(0,M^), with = 1, by the Helmholtz's decomposition Theorem (see 

Fujiwara-Morimoto [H] and Solonnikov [33]), there exist G E C°°(n) n Tyi'P'(O) and H e C°°{Tl)n 
LP'{n,M.^) such that 



9 
dG 

dv 



VG + H, divH = in 
g ■ 1/ on do,, 

C. 



Thus we have 



j V(div u)- Hdx = 



16 



so that 



j V(div u)-gdx = j V(div u) ■ (VG + H)dx = j V(div u) • VG < 

= /(V(divu)-^/).VGd. + ^//.VG<i. 

■ /" V X (curl u) ■VGdx + — ^ — /" f-VGdx 
J 2fi + X J 

f-VGdx, 



2n + X 
1 



+ A 



where we have used 



y V X (curl n) • VG = 0, 
since div(V x (curl n)) = in $7 and (curl n) x = on dil.. The above inequality implies 



V{div u)-gdx <||/||lp||VG||^,, <G||/|Up. 
Taking supremum over all such g^s, we obtain 

||V(div u)\\LP<C\\f\\L.. 

It is clear that this, with the help of (|3.3p and (|3.4p . implies (j3.2p . □ 
3.2 Existence 

In this subsection, we will first consider that (7 C is a bounded domain, and then employ the 
Galerkin's method to obtain a sequence of approximate solutions to (|l.ip - (|1.3p under ()1.6p and 
(jl.Sp or (jl.9p that enjoy a priori estimates obtained in §2, which will converge to a strong solution 
to (ll.ip - (ll.3p . The existence of strong solutions for the Cauchy problem on follows in a standard 
way from a priori estimates by the domain exhaustion technique, which will be sketched at the end 
of this subsection. 

To implement the Galerkin's method, we take the function space X to be either 

(i) for the Dirichlet boundary condition (jl.Sp . X := Hq n H'^{U,W^) and and its finite dimensional 
subspaces as 

X™ :=span{0\--- ,0"^}, m > 1, 

where {</>'"} C X is an orthonormal base of H^{Q,), formed by the set of eigenfunction of the Lame 
operator under the boundary condition n = on dQ; or 

(ii) for the Navier-slip boundary condition (ll.9p . 

X := {ue H^{n,R^) : u ■ u = {V x u) x u = on dn} , 

and its finite dimensional subspaces as 

X™ :=span{0\--- ,0"^}, m > 1, 

where {</>™} C X is an orthonormal base of H^{Q), formed by the set of eigenfunction of the Lame 
operator under the Navier-slip boundary condition u ■ u = (V x u) x u = on dQ,. By the W"^'^- 
estimate of Lame equation under (jl.Sp or (jl.9p (see Lemma 3.1), we see that {0™"} C W'^''p{VL) for 
any 1 < p < +oo. 



17 



Now we outline the Galerkin's scheme into several steps. 

Step 1 (modification of initial data). For (5 > 0, let Po = Po + dfj = do, and Uq G X be the unique 
solution of 



4 



C4 - V{P{pi)) - Ado ■ Vdo 
= 0; or Uq ■ v = {V X Uq) X u 



Po 9 in ^, 

on on. 



(3.5) 
(3.6) 



By the VK^'^-estimate of Lame equation, it is not hard to show that 



lim 

<5iO+ 



Uo - Uo 



X 



0. 



Step 2 (mth approximate solutions). Fix (5 > and 3 < g < 6. For m > 1 and some < T 
T{m) < +00 to be determined below, we let 



A;=l 



and look for the triple 

'p™ G C{[0,T];W'^^'i r\H^) 

m 

u"'{x,t) = E uf{t)(l)k{x) G C{[0,T];W^^^nH^) 
k=l 

G C{[0,T];H^n,S^)) 
solution of the following problem 
p^ + V ■ {p'^vJ^) = 0, 

(p"<, (t>k) + /i(V X n™, V(t>k) + {2p + A)(V • V<Pk) 
= -{p'^u' 
dY' + u^ ■ Vd" 
(p™, -u™, d"") 



■™ • Vn"^, (pk) - (V(P(p™)), </)fc) - (Ad™ • Vd'' 



{l<k<m) 



{u 



m ddr 



t=0 



Ad"" + |Vd""|2(i™, 
= (pi 4), 



(3.7) 



an X [o,T] 



0, or (n™ • z^, (V X u™) x 1^, ^) 



anx[o,T] 



0. 



The existence of a solution {p"^ , u'^ , d"^) to ()3.7p over Q x [0,T(m)] for some T(m) > can be 
obtained by the fixed point theorem, similar to that on the compressible Navier-Stokes equation 
by Padula [28] (see also [3]). Here we only sketch the argument. First, observe that for any given 
< T < +00 and u"^ G C([0, T]; W'^''^ n il^), it is standard to show that there exist 
(1) a solution p"^ G C{[0,T]; W^'" D H^) of ([321)1 along with p" 



Pi- 



(2) < < T, depending on n"' and ||do||i^3, and a solution d"* G C([0, t^], ^^([7^ 52^) ^ 



along with d™ 

It is well-known (cf. 



do and ^ 



0. 



er2x[o,t„ 

[3] or Lemma 2.5 in §2) that 



p™(x,t) > 5exp ^-^ ||Vu™||loo ds^ > 0, (x,t)GQT- 



(3.^ 
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The coefficients u^{t) can be determined by the following system of m first order ordinary differ- 
ential equations: 1 < k < m, 



(3.9) 



i=l 



where Ff. denotes the right hand side of p.7p 9. Since p"^ is strictly positive, the determinant of the 
m X m matrix {p^4>i-,4>k)i<_i k<m positive. Hence we can reduce (13. 9p into 



ut = Gu{uT,hT,t), bt = u^; uTiO) = 6^(0) = 0, 



(3.10) 



where Gk is a regular function of u^^h^. Therefore, by the standard existence theory of ordinary 
differential equations, we conclude that there exists a < Tm < tm and a solution vl^{t) to (|3.9p . 
which in turn implies the existence of solutions /o™, of (|3.7P i and (j3.7P ':i on the same time interval. 



Step 3 (a priori estimates). We will show that there exist < Tq < +oo and C > 0, depending 
only on the norms given by the regularity conditions on P and the initial data po,uo, and do, but 
independent of the parameters 6,m, and the size of the domain Q, such that there exists M ^ T 
so that for any m > 1, {(jf^ ,61^) satisfies: 



$'"(t) < exp 



Jo 



< t < To, 



(3.11) 



where ^"^{t) is defined by (|2.ip with {p,u,d) replaced by {p"^,u^,d"^) and 7W(po, '"o' ^o) is defined 
by (USD with {po,uo,do) replaced by (/o^, u^, d^). 

Since the argument to obtain ()3.1ip is almost identical to proof of Theorem 12. H we only birefly 
outline it here: 

First, it is easy to see (|3.7p 9 holds with 0^. replaced by n™. By multiplying ([321)3 by {Ad"^ + 
iVd^pd™) and integrating over Q and adding these two resulting equations, we can show that 
there is a M G such that the energy inequality (12. 5p holds with {p,u,d), M, and ^ replaced by 
{p"", -u™, d™), M, and 

Second, since (I3.7p 9 implies 



pm^jjn ^ V(P(/)™)) + Vd™ • Ad" 



(3.12) 



where ^m{u) = Yl'iLi{'^->'i^k)4'k ■ X — )• X"^ is the orthogonal projection map, we can check that the 
same argument as Lemma 12.31 vields that exists M £ so that 



Vn'"||^i < Af(«>"(t)), 0<t<Trr 



(3.13) 



Third, by differentiating (j3.12p w.r.t. t, multiplying the resulting equation with nj", integrating 
over ri, and repeating the proof of Lemma 12. 4^ we obtain that there exists M £ T such that for 
any m > 1, 



mi m 1 2 I 



Vm\2 ^ /~i 



Jn 



M{pi,u^,di)+ f M{¥%s))ds 
Jo 



Fourth, similar to the proof of Lemma 12.51 and Lemma 
that for all m > 1, 



(3.14) 

we have that there exists M £ F such 



C M{plu'^,dl)+ [ M{^^is))ds 
Jo 



(3.15) 
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and 



\p"^\\H^r^w^•1 < Cexp < C 



Jo 



(3.16) 



Fifth, by differentiating (j3.7p q w.r.t. x and mutiplying by Vd™ (and VAd™ respectively) and 
integrating over il., we can use the same argument as Lemma 12.71 and Lemma 12.81 to show that 
there exists M £ such that for all m > 1, 

||V2d™|||2+ r \\Vd'P\\l2ds<C[l+ f M{^"'{s))ds], (3.17) 
JO Jo 

||v3d'™||i2 + ^V'^^rili2t^s< (cA1(Po,<,4) + ^*Af($™(s))ds^ . (3.18) 

It is readily seen that combining all these estimates together yields p. lip with Tq replaced by Tm 
and Uq replaced by u™. 

Step 4 (convergence and solution). By the definition of Uq, A4 given by ()2.ip . and the condition 
pTTO]) . we have 

M{pl4,di) = 1 + 



and 



C 
< — 

- 6 



Uq -Uq 



0, as m — )■ oo. 



M{plu^,di)-M{pluldl] 

Thus there exists N = N{5) > such that 

M{pi,u^,di)<2 + \\g\\L2, ym>N. 

It follows from (|3.19p . (|3.1ip . and Gronwall's inequality (see, for example, [3] page 263 or [32 
Lemma 6) that there exists a small Tq > 0, independent of 5 and m, such that 



(3.19) 



sup <P"'{t) < C7exp(C||c/||i2), Vm > M. 

0<t<To 



(3.20) 



By virtue of p.20p . we obtain that for any m > M 



sup 

0<t<To 
(■To 



m m\\z I ||„m||2 i ||v7„."i||2 i ium||2 , ||v7j™i|2 



+ ^ ° {\\un\h„ + WVuTWh + W^'d^^Wh + W^'drWh) < Cexp{C\\g\\l,). (3.21) 

Based on the estimate p.2ip . we can deduce that after taking subsequences, there exists {p^ , ,d^) 
such that 

^ p^ weak* in L°^(0, Tq; W^''^ HH^), u"" ^ weak* in L°^(0, To;D^nD^), 

u"" weak in l2(0, Tq; D^'^), uf uf weak in L^{0, Tq; D^), 

V^u'P y^uf weak* in L"^ {0,T; L"^), 
^ / weak* in L°°(0, Tq; n D^) and ^^(0, Tq; L»^), 
dT 4 in L'^{0,T;H^) and weak* in L^{0,T;H^). 
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By the lower semicontinuity, (j3.2ip implies that for < t < Tq, {p^,u^,d^) satisfies 

:;5„<5ii2 



sup 

0<t<To 

rTo 







,<5||2 



74j<5||2 



2j<5||2 



(3.22) 



Furthermore, it is straightforward to check that {p ,u ,d ) is a strong solution in [0, Tq] of (II. ip - 
(jl.Sp under the initial condition {p^,u^,d^) = (pg,Ug,(io) ^^'^ the boundary condition (jl.Sp or 



t=o 



SJ. Since Tq > is independent of 6, {p^ , u^, d^) satisfies (j3.22p . p^ 



po in W^''inH\ uf, no 
in n D^, and dg = do, the same limiting process as above would imply that after taking a 
subsequence 510, {p^,u^,d^) converges (weakly in the corresponding spaces) to a strong solution 
{p,u,d) of (fm)- (fL3l) on O X [0,To] along with ([L6D and ([L8]) or (fLQl) . 

For the Cauchy problem on M'^, we proceed as follows. For R f oo, it is standard (cf. [Ml) that 
there exists dQ G //^(M^, S"^) such that = no outside Br for some constant no G 5"^ and 



lim 

i?,too 



0. 



Now we let Uq G Hq{Bji) n H'^{Bji) be the unique solution of 

- V(P(/?o)) - Ad^ • Vd^ = g on B^, 



as. 



0, 



(3.23) 



(3.24) 



where g G L^(]R^) is given by (jl.lOp . Extending Uq to by letting it be zero outside Br. Then 
it is not hard to show that for any compact subset K (ZM?, 



lim 

i?,too 



Vno 



m{K) 



0. 



(3.25) 



By the above existence, we know that there exists Tq > 0, independent of i?, and a strong 



solution {p 
condition: 



R „,R 



d^) of ([ri])-([L3]) on Br x [0,ro] of (iri])-([L3l), under the initial and boundary 



dd^ 



B.„.=o, = ^>">- ""'^ <""■ an. 



dBRx[0,To] 



0. 



{p^, u^, d^) 
Furthermore, {p^,u^,d^) satisfies the estimate: 

/TJio.Rw'^ _i_ 11^-^112 _L llY7-„-R||2 I IUR||2 1 ||Y7^-R||2 

/P Ut\\^2 + \\p |lH/i,9n//i + 11^^ llz/i + ll/fi + II vot 



(3.26) 



sup 

0<t<To 
To 



+ 



uYd^. + llVnflli. + ||V^d«||i. + llV^df Hi.) < Cexp(C7||<7||i2 



,R||2 



74jfl||2 



2 ,i?||2 



(3.27) 



with C > independent of R. It is readily seen that ()3.27p . ()3.23p . and ()3.25p imply that after taking 
a subsequence, we may assume that {p^,u^,d^) locally converges (weakly in the corresponding 
spaces) to a strong solution {p,u,d) of (jl.ip - (|1.3p on x [0, To] under the initial condition (|1.6p 
and the boundary condition (jl.7p . This completes the proof of Theorem 11.21 □ 
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(3.28) 



3.3 Uniqueness 

In this subsection, we will show the uniqueness of the local strong solutions obtained in Theorem 

Let {pi,Ui,di) {i = 1, 2) be two strong solutions on x (0, T] of ()l.ip - ()1.3p with (jl.6p and either 
()1.7p , or (jl.Sp , or ()1.9p . Set p = p2 — pi,u = U2 — ui,d = d2 — di. Then we have 

Pi + (ui • V)p + u ■ V/92 + pdivu2 + pidivu = 0, 
piut + piui • Vu + V {P{p2) - P{pi)) 

= Cu — 'p{u2t + U2 • VU2) — piu • Vu2 — Ad • Vd2 — Adi • Vd, 
dt- Ad = Vd- (yd2 + Vdi)di + \Vd2\^d - u ■ Vd2 - ui ■ Vd, 

with the initial condition: 

{p,u,d)\t=o =0, X G 0, 

and the boundary condition: 

Multiplying (I3.28P 9 by u, integrating over Q, and using integration by parts, we have 

^li J '^^'^'^ dx + J {{2p + A)|diviip + ^|V x up) dx 
= - J p{u2t + U2 ■VU2) -udx - J piu ■'Vu2 -udx + J {P{p2) - P{pi)) divudx 
+ j (Vd • VVd2 -u + Vd- Vd2 ■ Vu) dx - j Adi ■ Vd ■ udx. 
Observe that 

\P{P2) - P{pi)\ < Bp(\\p,\\L^ + ||p2||l-)|pI < C\p\. 



Hence, by Holder's inequality and Cauchy's inequality, we have 

pi\u\'^dx + j {{2p + A)|divnp + p\V x up) dx 



ld_ 
2di 



^ I|p||^|lk2t + n2 •Vu2||l6 11^11^6 + ||Vn2||L°° y^^il^l dx + \\p\\L2\\divu\\L2 

+ ||Vd||i2||V2d2||L3||u|U6 + ||Vd|U2||Vd2||L-||Vu|U2 + ||Adl||i3||Vd|U2||u| 

< ||p||^3||u2t + n2 • Vu2||L6||Vu||i2 + ||Vn2||vyi,9 j pi\u\'^dx 

+ ||p||L2||divu||i2 + ||Vd||i2||Vn||i2 

< e J |Vu|2 dx + C[||7J||% \\u2t + U2 ■ Vu2\\le 
+ \\Vu2\\w^„J pilnl'dx+WpWh + WVdWU. 



L6 
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Thus, by choosing e sufficiently small, we have 
d 



(3.29) 



^^jpi\u\^dx + j \Vu\^dx 
<C[||p||% \\u2t + U2 ■ Vnsllie + IIVusIIh^i,, j pi\u\^dx + W'pWh + ||Vd||i.j 
Multiplying (j3.28P i by 2/?, integrating over and using integration by parts, we have 
^ j IpI^ dx < j \'P^' ^P'2\ dx + j |^|^(|divtii| + |divu2|) dx + J \p pidivu| dx 

^llp||L2||V/92||L3||u||i6 + (lldiv^/illioo + lldivnallioo) j jpp + ||p||i2 ||divu||i2 
<llp|lL2||Vn||L2 + iWdivuiWwi.. + lldivTxalli^i,,) I \p\^dx (3.30) 
^llp||L2(||divM||i2 + ||V X u\\l2) + (||divni||vyi,9 + ||divu2||iyi,0 j \p\'^dx 
<e j \Vu\^ dx + CMlh + C(||divni||vi/i,, + lldivnall vi/i,q) / IpI dx, 
for any e > 0. Similarly, we have 

^ J \p\^ dx < J \p^ u-Vp2\dx + J \p\^ {\dwui\ + \divu2\) dx + J \p2 pidivu\ dx 



^l|p|l% l|Vp2||L2|kllL6 + (lldivuillLoo + ||divn2||Loo) / |p| 2 + ||p|| % ||div'a||^2 ||pi||^6 



<llnl|2 3 ||V'u||2,2 + (||divni||^yi,9 + ||divtt2||(4'i,g) / \p\^dx 



r\J 1 1 f 1 1 , 3. 



(3.31) 



^Ilpp3 I|Vu||l2 + (lldivnill^yi,, + ||divu2||(4'i,<j) / |p|2dx. 
]_ 

Multiplying (|3.31|) by ||p|| ^3 , and using Cauchy's inequality, we have 

3 < ||p||^3 ||Vllb2 + (lldiv^ill^,,, + ||divn2||H/M)||p||S 



<e y |Vupdx + C,||p||% +C(||divni||tyi,g + ||divti2||w'i,^ 



-||2 

3 • 



(3.32) 



Multiplying (j3.28P '^ by —Ad, integrating over fi, and using integration by parts and Cauchy's 
inequality, we have 

~ J \Vd\^dx + J |Ad|2dx<||Vd||i2||Ad||^2||V(i2 + Vdi||i^ + ||Ad||^2||Vd2|li6||d||i6 

+ ||Ad||i2||n||^6||Vd2||L3 + \\Ad\\L2\\ui\\L^\\Vdh2 

<||Vd|U2||Ad|U2||Vd2+Vdl|U2 + ||Ad||z.2||Vd2||^l||Vd||i2 

+ ||Ad||i2||Vnb2||Vd2|li2||Vd2|||6 + ||Ad||i2||Vni||^i||Vd||i2 

<||Ad||i2||Vdb2 + ||Ad||i2||V^Z||^2 

<i||Ad||i2+C||Vd||i2+C||Vn||i2. 
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This gives 

^ / \Vd\'^dx+ I \Ad\^dx<C\\Vd\\l2+C\\Vu\\l2. (3.33) 
dt Jn Jn 

Multiplying (j3.29p by 3C, putting the resulting inequality, (|3.30p and (|3.32p to (|3.33p . and taking 
e > small enough, we have 

^ (sCW^uWl, + IIpIIS + mh + llVdlli^) +C I \Vu\^dx 
^IIpII^I \\u2t + U2 ■ Vn2|||6 + ||Vti2||H/i.9 j pi\u\^dx + ||p||^2 + ||Vd||^2 

+(||divui||^yi,5 + ||divu2||vi/i.9) j \'p\^dx + ||p||% + (||divni||M/i,9 + ||divti2||w'i^ 
< {\\u2t + U2 ■ Vnsllie + IIVuiIIh/m + \\Vu2\\w^,, + l) 



(3.34) 



'iC\\^u\\l2 + \\-pfz + llpllia + l|Vd||i2 

Li 2 

By (I3.34p . Gronwall's inequality, and (po, ^^o, do) = 0, we have 

\\^/pl'u\\h + Wpf 3 + ||p||i2 + ||Vd|||2 + / / |Vnp dxds = 0. (3.35) 

Jo Jn 

This yields 

{p, u, Vd) = 0. (3.36) 
To see d = 0, observe that after substituting (j3.36p into ([328])3, we have 

dt = \Vd2fd, d\t=o = 0. 

This implies d = 0. This completes the proof. □ 

4 Proof of Theorem 11.31 

Let < < 00 be the maximum time for the existence of strong solution {p,u,d) to ()l.ip - ()1.3p . 
Namely, {p,u,d) is a strong solution to (|l.ip - (|1.3p in $7 x (0, T] for any < T < T^,, but not a 
strong solution in x (0, T*]. Suppose that (jl.lip were false, i.e. 

lim sup f||yO||Loo(o T-L°°) + / \\V d{t)\\la^ dt] = Mq < oo . (4.1) 
T/T. \ ' ' Jo ) 

The goal is to show that under the assumption (|4.ip . there is a bound C > depending only on 
Mq, po) 'UO) c?0) and such that 



sup 

0<t<T, 



maxdlpllpv-i,.- + 11/5*11^0 + (ll\/p^^tllL2 + l|Vn||^i) + {\dt\m + l|Vci||/^2; 

r=2,g 



< (4.2) 



and 



/ * (htll^i + ||7x||^2., + ||dt|||2 + ||Vd||^3) dt < C. (4.3) 

JO 

With (14. 2p and (14. 3p . we can then show without much difficulty that is not the maximum time, 
which is the desired contradiction. 

The proof is based on several Lemmas. 
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Lemma 4.1 Assume ( [^. we have 



[ * [ \V^d\'^dxdt < C. 
Jo Jn 



(4.4) 



Proof. To see (fO]) . observe that (jiJ]) implies /g^* ||V(i|||oc < Mq so that 



iVdl'^dxdt < Mo- sup l\Vd\^dx 



< M, 



P{p)\'^dxdt+ / (/OokoP + iVdoP) rfa; 



where we have used (j2.1ip in the last step. Applying (|2.1ip again, this then implies 

f ' f \Ad\'^dxdt = [ ' [ \Ad+\Vd\'^df dxdt+ [ * / \Vd\^dxdt 
Jo J Jo J Jo J 



< (1 + Mo) 



Jo 



P{p)\^dxdt+ j (pol^^ol^ + |Vdo|^) dx 



Since 



\P{p)\<Bp{\\pU^)\p\<C\p\, 
we have, by the conservation of mass and ()4.ip . 

/ / \P{p)\'^ dxdt < CT^, sup <C. 

Jo J 0<t<T, 

Thus the standard L^-estimate yields (j4.4p . 

Following the argument by [35], we let v = C^^V{P{p)) be the solution of the Lame system: 



□ 



\Cv = V{P{p)), 

\^\dQ ~ ^' or ^ ~^ \x\ — oo(when Q, - 
Then it follows from [35] Proposition 2.1 that 

W^vWl. < CWPip^L. < CBp{\\p\\lo^)\\p\\l. <C, l<q<6, 

where we have used (j4.ip and the conservation of mass in the last step. 
Denote w = u — v, then w satisfies 

pwt — Cw = pF — Vd • Ad, 
w\t=o = wo = uo- vo, 
w\q^ = or u) — )• 0, as |x| — )■ oo, 

where 

F = -u-Vu- C~^V{dt{P{p))) = -u-Vu + £-^Vdiv {P{p)u) - £-^V((P - P'{p)p)div u) . 
Then we have the following estimate. 



(4.5) 



(4.6) 



(4.7) 
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Lemma 4.2 Under the assumptions of Theorem \1.3[ if X < then {p,u,d) satisfies that for any 

o<t<n, 



8) 



/ {p\uf + \Vwf + \Vdf + \V^df)dx+ [ [ {{V'^^dl^ + \V^w\^ + \Vdt\^)dxds <C. (4. 
Jn Jo Jn 

Proof. The proof of this lemma is divided into five steps. 

Step 1. Est^mates of j dx. Multiplying gZJ, by integrating over ^, and usmg integration 

by parts and Cauchy's inequality, we have 



(;u|Vwp + + A)|div w\'^)dx + j p\wt\'^dx 

d f \ f ^ 

<\\^/pF\\l2+2—J{Vd0Vd--\Vd\%)■.Vwdx + C J \Vd\\Vdt\\Vw\dx = ^Ii. 



(4.9) 



For /i, we have 

h <\\Vpn ■ Vnlli. + W^C-'Vdiv (P(p)n)||i, + \\^C-'ViiP{p) - P\p)p)div u)\\l, 

' (4.10) 

i=i 

For 111, by Holder's inequality, (14. ip . Sobolev inequality, interpolation inequality, and (14. 6|) . we 
have 

1 1*6 
hi ^llpSulllsllV-ull^M < ||/55«|||6||V«||^2||V?X||26 



1 1 6 1 16 

<||p5n|||5||V'u||22||Vu;|||6 + H/JsuHis || Vn|||2 ||V?;|| ^ (4.11) 



4 / 6 
72„ 



<\\p'u\\is\\Vu\\l, \y\\V'w\\l, + ||Vu;||22 + Ij . 

Again by [35j Proposition 2.1, and (14. 7p . we have 

l|V't«||L2 < ||^/pu;t|U2 + ||Vpi^||i2 + ||Vd • Ad\\L2. (4.12) 

Substituting ()4.12p into (j4.1ip . and using Young's inequality, we obtain for any e > 

hi<e{\\^wt\\l, + \\^F\\l,) 

+C{\\p-su\\l4Vu\\l, + \\Vw\\l, + ||Vd||i^||Ad||i2 + 1). 

For Ii2 and I13, by [35] Proposition 2.1, (j4.ip . ()2.5p . and (jl.Sp . and Sobolev's inequality, we have 

^12 < ||^(p)n||i2 < \\pu\\l, < W^uWl, < C, (4.14) 

hs <IIV^IIi3||/:-'v((p(p) - p'(p)p)div u)\\l, 

<\\VC~'vmp) - P'ip)p)div u)\\l, < imp) - P'ip)p)Vu\\l, (4.15) 
<CBp{\\p\\lo.)\\p\\loo\\Vu\\12 < C\\Vu\\l„ 

where we have used the Sobolve inequality when i7 = M^, and both Sobolve and Poincare inequalities 
when ri is a bounded domain. 
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Putting (j4.13p . (j4.14p and (I4.15P into (|4.10p . and choosing e sufficiently small, we obtain 

h <\\\V~PWt\\l2 + C{\\p\u\\l4Vu\\l, + \\Vw\\l, + llVnIli. + \\Vd\\l^\\M\\l, + 1) 

<\\\VP^t\\l2 + C{\\plu\\U\Vu\\l, + \\Vv\\l, + llVnIli, + \\Vd\\l^\\/\d\\l, + 1) (4.16) 
<\\\^PMh + C(||p^u||i.||Vn||i. + llVnIli. + ||Vd||i.o||Ad||i. + 1), 



where we have used (j4.6p with q = 2. For L^, using Cauchy's inequality, we have 

h<- [ |Vdtpdx + C / |Vd|2|Vwp(ix 

^ (4.17) 
<^ j \Vdt\'^dx + C\\Vd\\l^ J \Vw\'^dx. 

Substituting (|i36|) and (jiTTl) into (jM]), we obtain 

^ J {p\Vwf + {fi + X)\diY wl"^) dx + ^ J p\wt\'^dx 

<2^ J{Vd0Vd-^\Vd\%) ■.Vwdx+^\\Vdt\\l2 (4-18) 
+c(\\Vd\\l.oi\\Vw\\l, + \\Ad\\l,) + ||p^n||i,||Vn||i, + ||Vn||i, + l' 



Step 2. EsUmates of J p\uf dx. Multiplying » by 5|n|3., integrating over n, and using inte- 
gration by parts and Cauchy's inequality, we have 



d 
It 



j p\u\^dx + j 5|np (/i|Vnp + (/u + A)|div np + 3^|V|n|p) dx 
= j 5P(p)div(|n|\)dx + j ^ {vd®Vd -]^\Vd\H^ div{\u\^u) - j 15(/u + A)(div u)\u\'^u ■ V\u\ 

<c(/p|.nv„i + / m'M) + /5(. + A,i«iv.v + / + A)i„nvNr. 

By Kato's inequality |Viip > |V|ti|p, we have 

r (15^ _ M/it^) ; |n|3|V|^z||2 > (15/i - ^^^) / |u|3|V^x|2, if ^ - ^(^^ < o 
\(15/i-^^^)/|n|3|V|^z||2 >0, if^-^^>0. 



Hence we obtain 



I J p\ufdx + 5 mm {p, (^4^ - ^i^^^ } | |u|3|Vt.|2dx 
<C{ p\uf\Vu\dx+ / \Vd\'^\uf\Vu\dx). 
Since A < we have 

Co := 5 min {p, (^Ap - ^il^±^^ | > Q. (4.20) 
27 



Thus by Cauchy's inequality, we have 

d_ 
di 



<C{j p\uf\Vu\dx + j \Vd\'^\uf\Vu\dx) 
<y j \u\^Wu\^dx + C j p^\u\^dx + j \Vd\'^\u\^dx 

Hence by Holder's inequality, Sobolev's inequality, the conservation of mass, (|4.ip and Young's 
inequality, we have 



^ j p\u\^dx + ^ j |np|Vn|2dx< j p^\u\'^ dx + j \Vd\ 

<(| p\ufdxy pUxy +\\V{\u\-i)\\l (I \Vdfdx 
<j j \u\^\Vu\'^dx + C 1 + J p\ufdx + ij \Vdfdx) 



u\ dx 



Thus by (I2.5P we have 



^ J p\ufdx + j J \uf\Vu\^dx< J p\ufdx+(^J \Vdfdx^ +1 



< 



p\ufdx + \\Vd\\loo\\Vd\\ls + I. 



gives 



(4,21) 



Step 3. Esi^ate. of j | V<i|= Differentiatmg (Q "th respect to we obtain 

Vdt - VAd + V{u ■ Vd) = V{\Vd\'^d). (4.22) 
Multiplying (|4.22p by SjVdpVd and integrating by parts over Q,, we have 

^ j \Vdfdx + 5 j \Ad\^\Vdfdx 

=5 J [V{\Vd\'^d) - V{u - Vd)] ■ \Vd\^Vddx - 5 j Ad ■ V{\Vd\^) ■ Vddx 

< J {\Vdf\V^d\ + \Vd\'^ + \Vu\\Vdf + \Vd\^\V^df) dx. 

This, combined with Cauchy's inequality and the fact 

\Vd\^ = -d- Ad (since \d\ = l), (4.23) 



^ J \Vdfdx + 5 J \Ad\'^\Vdfdx 

< I (iVdl^lV^dp + iVd^ + \Vu\\Vd\'\V^d\) dx (4.24) 

<||Vd||i.||v2d||i. + ||Vd||ioo||Vd||is + \\Vd\\U\Vuh4^'dh2 
<||Vd||i.(||V2ci||i, + llVnlliO + ||Vd||i.||Vd||i.. 
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By (USD and (023]), we have 

^ J \Vdfdx + 5 J \Adf\Vdfdx 

i\m\Ui\\v'd\\i, + \\vw\\i,) + iivdiii.. + iivdiii^iivdiiis. 



(4.25) 



<l 



Step 4. EsUmates of j dx. Multiplymg ^ by Vd„ integrating by parts over 1^, and 
using Cauchy's inequality, we have 



1 d 
2di 



\Ad\'^dx + J \Vdt\'^dx 
= I {V{\Vd\'^d) -V{u-Vd)) -Vdtdx 

(4.26) 

<e\\Vdt\\l2+C / (|Vd|'^ + |Vd|2|v2dp + |Vu|2|Vd|2 + |n|2|v2d|2)dx 
<e\\Vdt\\l2+c\pd\\l^{\\V^d\\l, + \\Vu\\l,)+ j \u\^\V^d\^dx], 
where we have used (|4.23p to estimate 

j \Vdf dx < II Vdllioo J I V^dl^ dx. (4.27) 



For the last term on the right hand side of (|4.26p . using Nirenberg's interpolation inequality and 
Cauchy's inequality, we have 



nl^lV^dl^dx 



5 

\U\ 2 



12 



" llVWso <e||V|n|t||i2+C||v2d||;" 



<e||V|tx|l||i,+C||Vd|||6||Vd||;|, 

<e||V|^|l||i2 + ellV^dlli, + C(||Vd||ia + WV^l^ + 1) 

<5e [ \u\^\Vu\^dx + e\\V^d\\l2+C{\\V^d\\l, + l). 



(4.28) 



By ()1.3p . ff^-estimate for elliptic equations, and (|4.27p . we have 
llV^dlli^ <||Vd,||i2 + ||V(^ • Vd)\\l, + ||V(|Vdpd)||2, 

<||Vd,||i2 + llVdIli.. (llVnIli. + llV^dll^,) + llVdIli. + / \u\'\V'd\'dx 



(4.29) 



<C[||Vd,||i2 + ||Vd||ioo (llVnlli^ + llV^dlliO + J \u\'\vMUx 

Substituting (|4.29p into (|4.28p . and choosing e sufficiently small, we have 

\u\'^\V^d\^dx <C j \u\^\Vu\'^ dx + llV^dll^j + ellVdill^a 
+ ||Vd||ioo(||V^||i2 + ||V2d||2^) + l 



(4.30) 
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Substituting (|4.30p into (j4.26p . using ()4.6p . and choosing e sufficiently small, we obtain 
d 



dt 
<C 

<c 



\M\'^dx + j \Vdt\^dx 



llVdIlioo (II V^dlli. + llVnIli,) + / Inl^l Vn|2 dx + \\V^d\\\, + 1 



llVdIlioo (||V2d||2, + II Vu;||i.) + II Vdlli^ + / |n|3| Vn|2 dx + \\V^d\\\, + 1 



(4.31) 



Step 5. Completion of proof of Lemma\4^ Adding (I4.21|) . (I4.18p . (I4.25P and M.Sip together, and 
choosing e sufficiently small, we obtain 



d 



j (p|n|5 + ^|Vu;|2 + (;u + A)|div w\'^ + \Vdt + lAdp) dx + ^ J p\wt\'^dx + \Vdt(' 



dx 



<2— {Vd(^Vd- -\Vd\%) -.Vwdx + C (llVnll^a + l) / plul^dxHVnll^^ + ||Vd|| 



dt _ 

+ llVdllioo (llVdIlis + llV^dlli. + ||Vu;||i.) + llVdIli. (||Vd||i5 + llV^dlli, + llVt^lliO 
+ llVdIlioo + llV^dlli. + 1^ . 

This, combined with Cauchy's inequality, implies 



j {p\uf + n\Vw\^ + (^ + A)|div + \Vdf + lAdp) dx + ^(^j p\wt\^ dx + 



d 
dt 

<2^ I {Vd^Vd- ^\Vd\%) -.Vwdx + C ||Vu||^2 + ||V(i||io 



\Vdt\'^dx 



+ + liv'rflli. + livdiiioo + 1) (iipi^iiis + livdilis + liv^diii. + iiv^iii.) + 1 

Integrating over (0, t), and using ()4.ip . ()2.5p . we have 

{p\uf + \Vw\'^ + \Vdf + \V^d\^)dx + J' J {p\wt\^ + \Vdt\'^) dxds 



<C 
where 



\Vd\'\Vw\ dx + / Kis) Wp^^uWI, + ||Vu;||i2 + llVdll^s + llV^dH^a ) ds + 1 



(4.32) 



K(.) = l|v^(s)|ii. + ||vMs)lli. + l|vd(.) 



+ 1. 



By (j4.32p and Young's inequality, we have 

{p\uf + \Vw\^ + \Vdf + \V^d\^)dx + J {p\wt\^ + \Vdtf) dxds 

<l J \Vw\^dx + c[J \Vd\Ux + K{s)(^\\p^u\\l, + \\Vw\\l2 + \\Vd\\ls + \\V^d\\l2) ds + 1 

<^( / I V«;|2 dx+ f \Vdf dx) + c[ r K{s) (\\p"^u\\l, + \\Vw\\l, + \\Vd\\l, + II V^dlli.) ds + 1 



Thus we obtain 



{p\uf + \Vw\^ + \Vdf + \V^d\^) dx + J (Hu-tp + \Vdt\^) dxds 



<C 



1 + / K{s) ( \\p-^u\\l, + llVti^lli. + ||Vd||^5 + llV^dlli. ) ds 
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(4.33) 



By dH]), ([23]) and ([13]), we know 

[ K{s)ds < C. (4.34) 
Jo 

By (I4.33p . (I4.34P and Gronwall's inequality, we obtain that for any < t < T^,, 

j [p\uf + |Vu;p + |Vci|^ + iV^dp) dx + j {p\wt\'^ + |Vdtp) dxds < C. 

This completes the proof of Lemma 14.21 □ 
Corollary 4.3 Under the same assumptions of Lemma\4.2[ we have that for any 2 < q < 6, 



sup (11^11^6 + ||Vn||i2 + llVdllL, + \\dt\\L2) + ||Vu||i2(o,T;L6) < C. (4.35) 

0<t<T. 



Proof. Combining (j4.6p with (j4.8p . we get 

l|Vn(t)||^2 < \\Vw{t)h2 + \\Vvit)U2 < C. (4.36) 
The upper bound of sup \\u\\i6 follows from (j4.36p and Sobolev's inequality. The bound of 

0<t<T* 

sup ||V(i||iQ follows from (14. Sp and interpolation inequality. For the last term of (I4.35p . by 

0<t<T. 

Sobolev's inequality, ()4.6p and ()4.8p . we have 

I|V'"IIl2{0,T;L6) ^l|V'w||z,2(o,T;Le) + II Vu||l2(o,T;L6) 

^l|V^U'||i2(o,T;L2) + ||Vu'||i2(o,T;L2) + 1 < C. 



By equation (jl.3p . (|4.8p and Holder's inequality, we have 

sup ||dt||L2 < sup {\\Ad\\L2 + \\Vd\\l, + \\u-Vd\\L2) 

o<t<r* o<t<T* 

< sup {\\uhe\\Vd\\Ls) + l<C. 
0<t<T* 

This completes the proof. □ 
Lemma 4.4 Under the same assumptions of Lemma \4-S\ {p,u,d) satisfies that for any <t < T^,, 

I {p\u{t)\'^ + |V(it|2)(t) dx+ f ( (iViip + \dtt\'^) dxds < C, (4.37) 

where f is the material derivative: 

f:=ft + u- V/. 

Proof. Step 1. Estimates of J p\u{t)\'^ dx. By the definition of material derivative, we can write 
(jl.2p as follows, 

pu + V(P(p)) =Cu-Vd- Ad. (4.38) 
Differentiating (I4.38P with respect to t and using (II. ip . we have 

M + />n • Vn + V(P(p)t) + (Vd • Ad)j 

r 1 (4-39) 

=Cu - C{u ■ Vti) + div Cu®u- V{P{p)) ®u - (Vd ■ Ad) u . 
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Multiplying (j4.39p by u, integrating by parts over Q and using the fact u = on dQ, we obtain 
/ '^1^1^'^^ + / (AilVup + (/i + A)|divnp) dx 
((P(p))tdivu + n(g) V(P(p)) : Vn) dx + fj, j {div (Au u) - A{u ■ Vu)) ■ iidx 

(4.40) 



+ (^ + A) J (^div{Vdivu0u) -Vdw{u-Vu)j ■udx + J {u {Ad ■ Vd)) : Viidx 

r ^ 

+ / (Vdt Vd + Vd Vdf - Vd • Vdth) : Vttdx = ^ Jj. 

1=1 

By equation (jl.ip and (|4.ip . we have 

Ji = y (^-div(P(/9)n)divit- (P'(p)p-P(p))divndivtt + ti® V(P(p)) : Vii^ 

= J {P{p)u ■ Vdivn + (P(p) - P'(p)p)divtidivn + P{p){Vuf : Vn - Pip)u ■ Vdivn^ dx 

= j {(Pip) - P'{p)p)divudWudx + P{p){Vuf : Vuj dx < \\Vu\\l2\\Vu\\l2 . 

By the product rule, we can see 

div {Au (g) n) — A(n • Vn) = Vfc(divuVfcti) — VkiV^u^Vju) — Vj(Vfcti-' V^n), 
so that by integration by parts, we have 

J2 = P j (Vfc(divnVfcu) - VkiVkU^Vju) - Vj{VkU^Vku)) -iidx < ||Vtt||i2 ||Vii|||4. 

Similarly, since 

div (Vdivn (g) u) - Vdiv {u ■ Vu) = V jU^ViU^) - VkiVju'ViU^) - \/i{Vku'VjU^), 
we have 

= + A) J {VkiVjU^Viu') - Vk(Vju'V,u^) - Vi(Vfcn'VjnJ')) tt'= dx < || Vn||i2 1| Vn||i4. 
By Holder's inequality, and Corollarv 14.31 we have 

J4 < ||Vn||i2||Ad||i6||Vd|U6||n||i6 < ||Vn||i2||Ad||i6, 



J5<l \Vu\\Vdt\\Vd\dx<\\Vuh4^dt\\L2\\Vd\\ 



Putting all these estimates into (j4.40p . using Young's inequality and Sobolev's inequality, and 
Lemma 4.2 and Corollary 4.2, we have 



p\u\'^dx+ / (;u|Vnp + (/i + A)|divttp) dx 



ld_ 

<||Vn||i2||Vtt||i2 + \\Vu\\14Vu\\l^ + ||Vn||i2||Ad||i6 + ||Vn||i2 ||Vdt||i2 ||Vd||Lo 
<f l|Vn||i2 + C {\\Vn\\l, + \\Vu\\l, + \\Ad\\l, + ||Vdi||i2||Vd||i.) 
<f l|Vn||i2 + C {\\Vn\\l + llV^dll^, + ||Vd,||i2 ||Vd||i. + l) 
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Thus we obtain 
d 



dt 



p\u\'^dx + fi / \Vu\'^dx<\\Vu\\l, + \\V^d\\l2 + \\Vdt\\l2\\Vd\\lo. + l. 



(4.41) 



By if'^-estimate of elliptic equations, Lemma 14.21 Corollary 14.31 and Nirenberg's interpolation 
inequality, we have 

\\^'dh2 <\\^dth2 + \\u • Vdh2 + |||Vn||V(i|||i2 + \\\Vd\\V'd\h2 + \\\Vdfh2 

<\Ndt\\L2 + \\u\\L4^d\\L3 + \\Vu\\L4Vdh^ + \\Vdh4^^dh3 + 1 



<l 



WVdth^ + \\Vu\\l4VurH, + ||V^(i||2,||V^d||^, + 1 

<^\\V^dh2 + C {\\Vdt\\L2 + \\V^u\\l2 + 1) . 

Thus we obtain 

WV^'dWL^ <\\Vdt\\L2 + \\V^u\\l2 + 1. 

By the definition ol w, we have 

Cw = pu + A(i • Vd. 

By ff^-estimate of the equation (I4.43|) . (14. ip . Corollary 14. 3^ Nirenberg's interpolation inequality, 
and ()4.42p . we obtain 



(4.42) 
(4.43) 



W'^Mh ^Wp^WI^ + W^d-vdWl, < \\ph\\l, + \\Ad\\Uvd\\% 



< 



\\p-2u\\l2 + \\Ad\\L2\\Ad\\Hi < \\p-^u\\h + \\VAd\\L2 + 1 



(4.44) 



<\\P~^u\\h + \\Vdt\\L2 + \\VML2 + l. 



By interpolation inequality. Corollary 14. 3^ (14.61) (for g = 6), (I4.44p . and Cauchy's inequality, we 
obtain 



l|Vn||i4 < 



< 



< 



< 
< 



\Vu 
\Vu 
\Vu 

\Vu 
\Vu 



{\\Vw\\l, + ||V^||i«) < iiv^xii^. iWvMh + 1) 

(llpUlli^ + ||Vdt||i2 + ||V2n||i2 + l) 

L4p'^u\\h + llV^^llie + INdtWl, + \\V'u\\l, + 1 
L4p'^u\\l, + \\Vdt\\l2 + \\V'u\\l, + l. 



(4.45) 



Putting ^^J5\i and (j02|) into (031]), we have 

d_ 

di 



(4.47) 



I p\uf dx + pj I Vn|2 dx < \\Vu\\l, llpkllia + llVdt lli^dl Vdllioo + 1) + HV^Hi^ + 1. (4.46) 

Ste, 2. EsUmates of J | V..|^ dx. Differentiating » with respect to t, we have 

dtt - Adt = dt (I Vdpd -u-Vd) . 
Multiplying ()4.47p by du, integrating by parts over and using ^ = 0, we obtain 

^ ^ ^ iVdfp dx + j Idttp dx = j dt (iVdpd -u-Vd) dtt dx 

(iVdpldtl + \Vd\\Vdt\) \dtt\ dx + j {\ut\\Vd\ + \u\\Vdt\) \dtt\ dx = Ki + K2. 



2dt 



< 



(4.48) 



33 



By Holder's inequality, Sobolev's inequality, Corollary 14.31 and Young's inequality, we have 

l^il <\\du\\L4dt\\L4^d\\l, + ||d«||i2||Vdt||i2||Vd|Uoc 
<\\dtt\\L<\\^dt\\L2 + 1) + \\dtt\\L4'^dth2\\Vd\\L^ 

<l\\du\\h + C (llVdilli. + ||Vd||ioo||Vd,||i. + 1) . 

By the definition of u, Holder's inequality, Sobolev's inequality. Corollary 14.31 and Young's inequal- 
ity, we have 

\K2\ < J [{\u\ + \u\\Vu\)\Vd\ + \u\\Vdt\] \du\ dx 

<\\dtt\\LA\u\\LA\^ d\\L3 + ||dtt||L2 11^11^6 llV-ull^ellVfillie + ||rftt||L2 ||n||i6 1| Vdj H^s 



<\\dtt\\LAWu\\L^ + \\dtt\\L^{\\V\\\L2 + 1) + \\dtt\\LAWdt\\L-^ 

<\\\du\\h + C {\\Vu\\l2 + + l|Vd||i3 + 1) 



Putting these two estimates into (j4.48p . using Nirenberg's interpolation inequality, and Young's 
inequality, we have 



\^J\Vdt?dx + \j\du?dx 



<livu|ii2 + iiv\|ii2 + (1 + \\yd\\i^)\\vd,\\i, + iivdilia + 1 

<|IV7„M|2 , IIV72„,||2 



(4.49) 



llV^^lli^ + \\V'u\\i, + (1 + \\Vd\\U)\\Vdt\\i, + \\Vdt\\L2\\V'dt\\L2 + 1 

<\\\V^dt\\l, + C {\\Vu\\l, + + (1 + l|Vd||i..)||Vdi||i2 + 1) . 

By ff^-estimate of the equation (|4.47|) and estimates similar to Ki and i^2) we obtain 

\W^dt\\L2 <\\dtt\\L^ + \\dt{u ■ Vd)\\L2 + \\dt{\Vd\^d)\\L2 

<\\dtt\\L^ + \\u ■ Vd\\L2 + \\{u ■ Vn) • Vd\\L2 + \\u\\LA\^dt\\L^ 

+ Mt||L«l|Vd||ie + ||Vdi|U3||Vdb6 
<\\dtt\\L^ + ||^||L«l|Vd||z.3 + h||i6||Vn||i6||Vd||i6 + ||Vdi|||,||Vdt|||e + \\Vdt\\L2 + 1 

<\\\y^dt\\L2 + C {\\du\\L2 + ||Vn||i2 + + l|Vdt||L2 + 1) . 

Thus 

l|V'dt||i2 <\\dtt\\L^ + ||Vn||i2 + + l|Vdt||i2 + 1. (4.50) 

Substituting this inequality into (|4.49p . we obtain 

j \Vdt\''dx + I |<i«|2dx<||Vu||i, + ||V2n||i2 + (l + ||Vd||ioo)||Vdi||i2 + l. (4.51) 



d_ 

di 



Combining ()4.46p and (j4.5ip . and applying Gronwall's inequality, we establish the conclusions 
of Lemma 14. 4i □ 

By the equation (j4.43p and Lemma 14.41 we obtain the following Corollary. 
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Corollary 4.5 Under the same assumptions of Lemma\4-S\ we have that for q E (3,6], 



sup [\W^d\\L2 + \\Vd\\L^) + ||Vu;||i2(o,T.;L-) + \\^ M\mo,T.;L'>) < C (4.52) 

o<t<n 

Proof. By '^-estimate of elliptic equations, (|1.3p . Lemma 14.41 Corollary 14.31 and Nirenberg's 
interpolation inequality, we have 

\N''d\\L2 <\Ndt\\L^ + \\u-Vd\\L2 + |||Vn||Vd|||i2 + \\\Vd\\V^d\\\L2 + \\\Vdf\\L^ 

<||n||i6||Vd||i3 + \\Vu\\L2\\Vd\\L^ + ||Vd||i6||V2d||i3 + 1 

<l|vd|lil|vd|||, + wv'diiwvmi, + 1 < ^WV^Il^ + C. 

Hence 

sup \\V^d\\L2 < C. 
o<t<r» 

By Sobolev's inequality, this yields 

sup llVdllioo < C. 

0<t<T, 

For simplicity, we only consider the case q = 6. By VF^''?-estimate of the equation (j4.43p . (j4.ip . and 
Sobolev's inequality, we obtain 

W'^^MIl^ <\\pu\\l<^ + II Ad • Vd\\L6 < \\u\\l6 + \\Ad\\Hi II V^IIl- < \\Vu\\l2 + 1. 
Therefore, by (I4.37p . we have 

l|V2w;||i2(o,T.;L6) < / * (llVnIlia + l) ds < C. 
Jo 



□ 



sup \\Vp\\l,^l2<C. (4.53) 

0<t<T, 



Following the same argument of [35] Section 5, we have 
Lemma 4.6 Under the same assumptions of Lemma\4.S\ we have that for q £ (3,6], 



Corollary 4.7 Under the same assumptions of Lemma\4.2[ we have for q G (3,6] 



sup ||V^n||i2 + ||u||^2(ot,-d2.'3) < C"- (4.54) 

0<t<T, 

Proof. By Proposition 2.1 in [35], (|4.38p . (|4.ip and Lemma 14.61 we obtain that for ri = 2 or q, 

llV^ullin <\\pu\\Ln + ||V(P(p))||l.i + IlVd • AdWiri ^ , 

(4.55) 

<||pn||Ln + II Vd • Ad||L'-i + llVpllin . 
When ri = 2, (14. ip . (I4.55p . Lemma [42| Lemma [44], and Corollary 14.51 imply 

I|V\||l2 < ||pk||i2 + ||Vd||L^||Ad||i2 +1<C. 

When ri = q, for simplicity, we only consider the case q = 6. By (14. ip . (I4.55p . Lemma l4.2[ Lemma 
14.41 Corollarv 14.51 and Sobolev's inequality, we have 

I|V^'"IIl2(0,T.;L6) ^l|p||L°°(0,T.;L°°)ll^llL2(o,T.;Le) + SUp || V(i||Loo || Ad||i2(o ^^^^6) + 1 

o<t<n 

^I|V'u||l2(o,t.;L2) + ||A(i||i2(o,T.;Hi) + 1 < C. 

This completes the proof. □ 
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Corollary 4.8 Under the same assumptions of Lemma\4-S\ we have that for ri = 2 or q, 



sup {p\ut\^ + \ptr)dx+ {\Vut\^ + \V^dt\^ + \V^d\^) dxds<C. (4.56) 

o<i<T, Jn Jo Jn 



Proof. It follows from (14. ip . Lemma 14.41 Sobolev's inequality, (14.35p . and Corollary 14.71 that 
p\u,\^d.<lp\u\^dx + lp\u.Vu\Ux 

<\\p\\l^\\u\\loo j \Vu\''dx + l <\\Vu\\h^ + 1<C. 

By (fTT]) . Sobolev's inequality, (j05]) . Lemma SSI and Corollary we get 

ll/OtlU'-i ^ll/odivn||Lri + ||u • Vp||l'-i < ||/9||lcx) ||divn||Lri + || V/oHLn 

<||/>||L-||divn||^i + ||V^x||^i||Vp|U'-i <C. 

By Lemma 14.41 interpolation inequality, Sobolev's inequality, (j4.35p . and Corollarv 14.71 we have 



Jn 



|Vnt| dxds 



\Vu\ dxds + 



Jn 



|V(u • Vti)r(ix(is 



Jn 



< I I \Vu\^dxds+ I I \u-V^u\^ dxds + l 
Jo Jn Jo Jn 

< I \\Vu\\L2\\Vu\\]j^ds + I \\u\\\oo I |V\pdxcis + l 
Jo Jo Jn 



^ / l|Vn|||^i / \V'u\' dxds + l<C. 
Jo Jn 



By ()4.50p , Lemma 14.41 aiid Corollary 14.71 we get 

r-T* 



\VMt\ dxds < C. 

'0 Jn 

By ff^-estimate of the equation (|1.3p , we have 

\N'd\\l, < WdtWl. + \\n- Vd\\l, + \\\Vd\'d\\l, = J2 L^. 
For Li, ()4.35p and Lemma 14.41 implv 

Li<\\V^dt\\l2 + l. 



(4.57) 



i=l 



(4.58) 



(4.59) 



For L2, Holder's inequality, Sobolev's inequality, (|2.5p . (|4.8p . (|4.35p . Corollarv 14.51 and Corollary 
14.71 we have 



L2 



nldVdl + |V2(i| + |V='(i|) 

+ I|V^* + ||V 
+ \\Vu\\UNmir<C. 



|Vn|(|V(i| + iV^dl) + \V^u\\Vd\ 



L2 



L2 



<h||io. ||Vd"2 , 11V72^„2 , „V73^„2 \ , „w^„2 /„w.,„2 , „w2.,„2 



+ ||Vd||i^ ||Vn||i. + ||V^n||i. 



(4.60) 
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Similarly, for L3, we have 

L3 <ll|vdp||2, + \\\s7d\\v'd\\\l, + \\\v'd\YL^ + iiivdiiv^diiii, 

<\\Vd\\U\Vd\\j,, + \\V'd\\%,<C. 
Substituting (039]) -(jMU) into K5E\i . we have 

\\^'d\\h<\\^'dt\\l, + l. (4.62) 
Integrating (j4.62p over (0,t), and using (|4.57p . we establish Corollary 14. 8[ □ 

Proof of Theorem [TSl 

By the above estimates, we know that both (|4.2p and (|4.3p are valid. Hence T^, is not the 
maximum time for the strong solution {p,u,d). This contradicts the definition of T^. The proof of 
Theorem 11.31 is complete. □ 
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